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Charge reversal of the planar electrical double layer is studied by means of a well known integral 
equations theory. By a numerical analysis, a diagram is constructed with the onset points of charge 
reversal in the space of the fundamental variables of the system. Within this diagram two regimes 
of charge reversal are identified, referred to as oscillatory and non oscillatory. We found that these 
two regimes can be distinguished through a simple formula. Furthermore, a symmetry between 
electrostatic and size correlations in charge reversal is exhibited. The agreement of our results with 
other theories and molecular simulations data is discussed. 

PACS numbers: 



I. INTRODUCTION 

Charged particles are naturally adsorbed onto an op- 
positely charged surface, however, temperature prevents 
them to completely condensate producing a diffuse ionic 
concentration profile known as the electrical double layer 
(EDL). Intuitively one might expect that the adsorbed 
counterions in the EDL are just the necessary to com- 
pensate the surface's charge, however, for certain condi- 
tions, they overcompensate it producing a surface charge 
reversal. Hence, an inversion of the local electric field 
takes place next to the first counterions layer and, con- 
sequently, coions are attracted to form a second layer 
referred to as charge inversion of the EDL. 

Charge reversal (CR) and charge inversion (CI) have 
motivated a large number of studies in the past [l], 0, 
S, S H S 0, 1, H, 01 [U Gl M, M El. These 
effects have been observed in the formation of self- 
assembled polyelectrolyte layers on a charged substrate 
[ill, self- assembled DNA- lipid membrane complexes 17 1 
and anomalous macroions adsorption on Lagmuir films 
[Tc| . Furthermore, CR has been associated with re- 
verse of the electrophoretic mobility of charged colloids 
[H[2(J[HL[i3 and polyelectrolytes [23j| as well as with at- 
traction between like ch arge d macroparticles 0, HH, HH 
and DNA condensation [27j |. 

CR and CI manifest as alternated layers of coions 
and counterions in the neighborhood of a charged in- 
terface, i.e., imply spatial oscillations of the inhomogc- 
neous charge distribution profiles. For bulk electrolytes, 
the oscillatory behavior of the pair correlation functions 
was predicted by Stillinger and Lovett [28| as a conse- 
quence of their second moment condition. The oscilla- 
tions of the inhomogeneous charge density profiles have 
been observed in studies where the restricted primitive 
model (RPM) electrolyte is used to model the diffuse 
EDL 0, [H, [H U m, H, [H, H|. For the EDL pro- 
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duced by a RPM electrolyte at a planar interface, two 
regimes of CR have been distinguished [HI, HH|, i.e., (i) 
the oscillatory and (ii) non oscillatory. These regimes are 
named according to the behavior of the electrolyte pair 
correlation functions in bulk. In the former regime, since 
the bulk pair correlation functions oscillate, any distur- 
bance caused by an external electric field is propagated in 
the same way, therefore, CR occurs for any nonzero sur- 
face charge density. In the second regime, a sufficiently 
high surface charge density is necessary to produce CR. 
In both regimes, short range and long range correlations 
correlations play an important role in the surface CR 

am. 

For some cases, the onset points of oscillations of the 
bulk ionic pair correlation functions [3(| and EDL-CR 
[1, @, have been computed, however, a complete pic- 
ture of both phenomena and their interrelationship is still 
missing. To fill this void, we construct a CR diagram for 
the planar EDL in terms of its fundamental variables, 
i.e., the particles volume fraction (77), the ion- ion (£) and 
ion-wall (7) Coulomb coupling. We also analyze the de- 
pendence of CR on each variable. Hence, the reminder 
of the paper is organized as follows: in section [IT] we de- 
scribe the model and theory. Section IIIII contains the re- 
sults and their discussion, and, finally, some conclusions 
are given in the closing part. 



II. THEORY 

Particles and fields are equivalent since both are de- 
fined thorough their interaction with other particles and 
fields. This simple fact has been used in physics of in- 
homogeneous fluids to derive integral equations theories 
[38[, i.e., in a nonuniform fluid the external field is taken 
as one more species of the fluid. Hence, one can use stan- 
dard integral equations theories for uniform fluids based 
on the Ornstein-Zernike equation. 

In the Ornstein-Zernike equation for inhomogeneous 
fluids, the total correlation function between an external 
field and a particle of the fluid, h a i (r2i), is decomposed 



Typeset by REVTeX 



2 



r 




FIG. 1: Set up of the system under study. 



as the sum of direct and indirect correlations, i. e., 

n „ 

h a i(v2\) — Cgj (r2i) + ^""^ / p m h am (r23)c m i(r 3 i)dr 3 (1) 

where c Q i(r 2 i) represents the direct correlation function 
between the external field and a particle of species i lo- 
cated at Y21 and c m j(r3i) is the direct correlation func- 
tion between any pair of particles, of species m and i, 
in the fluid. In this notation the subindex a refers to 
the external field. The second term in Eq. JT]) represents 
the indirect correlation between the external field with a 
particle of species i, carried out through all the particles 
of species m in the fluid. The total correlation func- 
tions are related to the particles distribution functions, 
ffmfoi), by h m (r 2 i) = 5m( r 2i) - 1. 

The Ornstein-Zernike equation requires extra rela- 
tionships, referred to as closures, between h a i{v2i) and 
c-ai ( r 2i) to be solved. Here, we apply the hypernet- 
ted chain/mean spherical approximation (HNC/MSA), 
which implies that c Q ;(r2i) is approximated by the hy- 
pcrnctted chain (HNC) closure, whereas c m j(r3i) is taken 
through the mean spherical approximation (MSA). This 
approach has been applied for modelling inhomogeneous 
and confined electrolytes in many geometries and, in gen- 
eral, it has shown good agreement with computer simu- 
lations USEE!- 

A. The surface and fluid models 

We consider an electrolyte solution contiguous to a pla- 
nar wall which has a uniform surface charge density <tq 
and dielectric constant e w (See Fig. [T]). In the restricted 
primitive model (RPM) electrolyte, the ions are consid- 
ered as hard spheres of diameter a, with a centered point 
charge qi = ezi(e is the proton's charge and z\ is the ionic 
valence), and at concentration pi. The ions are embed- 
ded in a structureless solvent, considered only through a 



uniform medium of dielectric constant e and, for simplic- 
ity, e w = e. Within this model, a high particles volume 
fraction (large a and/or high pi) could be a simplistic way 
of taking into account solvent effects, such as solvation 
or excluded volume. In this regard, some recent works 
have considered the solvent effect through more refined 
models (UIH- 

In planar geometry, g a i(r2i) depends only of the per- 
pendicular distance to the surface, x, therefore here- 
inafter we will simply write giix) = g a i(?2i) and h m (y) = 
h a m(j23), where the subindex a has been omitted. The 
charge on the wall is compensated by an excess of charge 
in the fluid given by the charge density profile, p e i(x), 
such that 



- CT = / Pei(x)dx (2) 

J a/2 

The charge density profile gives the structure of the elec- 
tric double layer (EDL) and it is expressed in terms of 
the ions distribution functions, gi{x), as 

2 

p e l(x) = ^2 ziepigiix) (3) 

i=l 

The surface-ion direct interaction potential is written 
as u a i{x) = Uai(x) + where the first term is the 

ion-surface hard-core interaction potential which consid- 
ers that ions can not penetrate or deform the surface. 
The second term is the electrostatic contribution, given 
by 

27T 

u ai( x ) = —qiVQ(x ~ Xoo) (4) 

being Xoo a reference point. 

The i species concentration profile, gi(x), is related to 
the ion- wall potential of mean force, Wi(x), by 

gi(x) = exp{-/3wi(x)} (5) 

from which the ion- wall effective mean force, f x i{%), is 
obtained by 

dwjjx) d\ng t (x) 
Jxi(x) = 7^— = k B T — (6) 

From this relationship one can infer that a negative (pos- 
itive) slope of gi(x) implies an attractive (repulsive) be- 
havior of fxi(x). 

The computation of gi{x) requires to solve Eq. |T]). 
As we pointed out above, it is solved by approximat- 
ing Ci(x) = c a i(r2i) by the hypernetted chain closure in 
Eq. {T}, i.e., 

Ci(x) = —]ngi(x) — /3ui(x) + hi(x) (7) 
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Thus, by combining Eqs. J7]) and |T]) we get 

gi{x) = cxp { -pu ai (x) + p, 

m— 1 



h m (y)cl SA (s) dv 3 



(8) 

In addition, in Eq. ([8]), we used the bulk MSA direct 
correlation functions between two ions of species i and 
m, i.e., Cmi' A (s), with s=|r3i| their relative separation 
distance. As usual, /3 = 1/(/cbT), where fcs is the Boltz- 
mann constant and T is the absolute temperature. The 
explicit form of c^f A (s) for the RPM electrolyte can be 
written as 



MSA i 



(s) = c hs (s)+q m q iC ST (s) 



Pq m qi 



(9) 



where c hs (s) and c sr (s) are short range functions, being 
equal to zero for s > a. By integrating Eq. ||SJ) in cylin- 
drical coordinates, after some algebra it is written as 



9i{x) 



exp ■ 



47T/3 



Zieagx + 2ir 



+ J(x) + 2-kz 
2iTf3e 2 z. 



la/2 

hd{y)L(x,y)dy 

a/2 

hd{y)[x + y 



a/2 




with i = 1,2 and defining h s (y) 
h d{y) = ELi z mPmh m {y) and 



T, m =lPm h m(y), 



K(x,y) 
L{x,y) = / 

J \x 

J(x) 



c s (s)sds 



c d r (s)sds 



a/2 



K(x 7 y)dy 



In Eq. (|13p size correlations are included through the 
first two terms, meanwhile, the last term includes elec- 
trostatic short range correlations. Thus, in the case of 
qi = 0, ip(x) = and only the size correlations terms 
are included in Eq. (fTTj) . hence, the model reduces to the 
inhomogeneous hard-sphere fluid. On the other hand, in 
the point ions limit (a = 0), wf(x) = and Eq. (jTT|) 
becomes 



i(x) 



exp{-0qiijj(x)} 
2tt 



cxp 



4-7T 



-(Jo(3q l x 



p c \(y)[x + y + \x-y\]dy } (14) 



a/2 



which is the integral version of the Poisson-Boltzmann 
(PB) equation. In Eq. (fT4)l a closest approach distance 
of the ions to the wall, a/2, is considered. PB equation for 
size-symmetric electrolytes does not predict CR Q since 
both, size and electrostatic short range correlations are 
neglected. However, it should be mentioned that models 
based on the PB equation may predict charge reversal 
when size correlations are considered artificially, i.e., by 
assigning to cations and anions unequal closest approach 
distances to the wall [43[ . 

Eqs. pT|) and (fTl)) depend on several parameters, i.e., 
on pi, a, T, e, Zie, <7q. To simplify the analysis below 
we introduce the three fundamental parameters that de- 
scribe the planar EDL. 



B. The fundamental parameters 

The thcrmodynamical properties of the symmetric 
RPM electrolyte can be described exclusively in terms 
of the following parameters [28[ 



ae 



(15) 



where px = X)j=i Pi ano - ^ has been used that hj(x) = 
-1 for x < a/2. The HNC/MSA integral equations, 
Eq. (I10p . can be conveniently written as 

g.i(x) = exp{-/3wi(a;)} = exp {-f3qi^{x) + wf (x)} 

(11) 

where ip( x ) is th e mean electrostatic potential, expressed 
as, 



4tt 2tt 

ip{x) = a x 

e e 



a/2 



Pe\{y)[x + y+ \x-y\]dx (12) 



and wY(x) is a short range potential, given by, 



n 



7T 

6 



Pra" 



(16) 



^ accounts for strength of the ion-ion electrostatic inter- 
action and is referred to as the ion-ion Coulomb coupling. 
The second parameter, 77, is the particles volume fraction. 

The ion-wall direct electrostatic interaction at the 
wall's surface is given by 



Ui 



a\ 2nqiao 



a 

X °° ~ 2 



thus, we define a third parameter 

27T/3 



li 



-qi<j a 



(17) 



(18) 



wf(x) = -2vr / h s (y)K(x,y)dy-J(x)-2m 

J a/2 



^hjch^^iantifics the ion- wall electrostatic interaction. 
a / 2 We 'also define 7 = |7j|. It is an easy task to show that 
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total potential energy and the partition function of the 
system (charged wall + electrolyte) can be written in 
terms of these dimcnsionless variables. 

In terms of the 7$ parameter, Eq.Q can be written as 



(3u m {X) = ^(X-Xoo) 



(19) 



with X = 2x/a the dimcnsionless distance. Hence, 
Eq. (fT0|) is written as 



9i (X) = exp| 74 X + J h as {Y)K{X,Y)dY 

h ad (Y)L(X,Y)dY (20) 

/OO 
h ad (Y)[X + Y + \X-Y\]dY 

where Zj = Qj/q and we have defined 

n 

h as (X) = ^2 Vmham(X) 
m—1 

n 

h a d(X) = ^2 z m -q m h am (X) 



4tt 



E(x) = —[a + a'(x)] 



(21) 



being er'(x) the surface charge density accumulated be- 
tween two planes, parallel to the surface at a/2 and x, 
respectively, given by 



a'(x) 



a/2 



Pv\{y)dy 



(22) 



The electroneutrality condition implies that 
lim <r'(x) = — f7o, such that the electric field pro- 

x — * 00 

duced by the plate is completely screened at infinity. 
In order to evaluate charge reversal (CR) we define the 
dimcnsionless function 



a(x) 



cr + a'(x) 



(23) 



such that CR occurs in an interval [0:1,352] where a'(x) 
overcompensates a Q , or cquivalcntly, a(x) > 0. Conve- 
niently we also define a* as the absolute maximum of 
<r(x), i.e., 



The expressions for K(X,Y), L(X,Y) and J(X) are 
given below. ForX-2<F<X + 2, 

K(X,Y) = J Cl [(X-r) 2 -4] + ^ C2 [|X-r| 3 -8] 



160 



??ci[|Jf-F| 5 -32] 



L(X,Y) = 3[2-\X-Y\]-l-^— [A-(X-Yf] 



Ta 



2(1 + Ta) 
\X-Yf] 



4 \1 + IX 

otherwise K(X,Y) = L(X,Y) = 0. For 1 < X < 3 

J(X) = i ?7 c 1 [p 3 -12 J 5+16] + |?/ 2 c 2 [p 4 -32p + 48] 
4 o 

+ ^ 2 ci{P 6 - 192?3 + 320] 

otherwise J(X) = 0, with p = X - 1. Eq. J20]) is solved 
using a finite element technique. Details of the numerical 
method are given elsewhere [Hj]. We wish to point out 
that the reduced concentration profiles, gi{X), depend 
only on the fundamental parameters £, 77 and 7, for any 
combination of pi, a, T, e, ^e, a . That is, two systems 
are equivalent provided that 77, £ and 7 are the same for 
both systems. 



III. RESULTS AND DISCUSSION 

The local electric field, E(x), is useful to identify and 
quantify charge reversal. By Gauss' law it is given by 



a* = max{cr(x) : x £ [a/2, 00)} (24) 
which, from Eq. |22|) . is written as 

1 f 1 "- 

a = -1 / p e i(x)dx (25) 

^0 Ja/2 

being x rna x the location of the absolute maximum of 
a(x). For the symmetric RPM electrolyte, which is the 
case considered here, z + = —z^~ z and p+= p—= p s - 
Furthermore, it can be readily shown that the extreme 
values of <j(x) are located at a set of points, Xi, where 
Pei(xi) = J2i=i z i e Pi9i(%i) = 0, implying g+(xi) = 
g—(xi). It will be shown that the absolute maximum 
of a{x), a*, is the closest one to the wall. Thus, the 
plane at x — x max defines the boundary of the adsorbed 
countcrions layer and a* quantifies the rate of charge 
overcompensation. 

We constructed a charge reversal diaaram\4^ in the 
77-^-7 space in the following way: For a given constant 
value of 7 and starting from values of £ and r] close to 
zero, both parameters were progressively raised until the 
condition a* > was reached (for numerical purposes we 
set a* > 10~ 5 ). Hence, by repeating this procedure for 
different values of 7 the full space was analyzed. The on- 
set of points of CR (//*,£') define a curve referred to as the 
charge reversal curve. Thus, the main result of this work 
is presented in Fig. [21 the charge reversal diagram for 
the planar EDL. For clarity only three 7-dependent CR 
curves are plotted. For each value of 7, below the curve 
charge reversal does not occur whereas it does above it. 
It should be mentioned that a* continues increasing for 
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Non-oscillatory regime 



0.02 ' 0.04 ' 0.06 ' 0.08 ' 0.1 



FIG. 2: Charge reversal (CR) diagram in the £-77 plane with 
three CR curves: 7 < 1 (solid line), 7 = 5 (dashed line), and 
7 = 8 (dotted line). The systems described in Table I are 
plotted within the diagram: Black symbols are onset points 
of CR, whereas white and light symbols mean the absence and 
occurrence of CR, respectively. The inset is an amplification 
for high values of £. For the onset points of CR from molecular 
simulation (b,c,f) the numerical error is approximately the 
symbol size. 
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TABLE I: Numerical values (77 and £ ) for the points (a-h) 
in Fig. [2] with their corresponding values of ionic diame- 
ter, a, electrolyte concentration, p s , and valence, z. In all 
cases Ib = ~ 7.14A. Method refers to dressed ion the- 
ory (DIT) , molecular dynamics (MD) and Monte Carlo (MC) 
simulations. 



I (£77 > 0.0662) and regime II (£77 < 0.0662). In regime 
I, CR occurs for any value of 7 different from zero, that 
is, CR is a property of the fluid due to the fact that the 
bulk ion-ion pair correlation function is oscillatory. This 
implies that any perturbation in the fluid (independently 
of the external field strength) is propagated in an oscilla- 
tory manner, producing alternated layers of counterfoils 
and coions next to the wall, hence, CR and CI. Although 
not shown, we found that this is valid for different geome- 
tries of the external field. In regime II, the ion-ion pair 
correlation function is monotonic, hence, the occurrence 
of CR depends on the external field. According to the 
behavior of the pair correlation function in bulk, regimes 
I and II are referred to as oscillatory and non-oscillatory, 
respectively. 

The onset points of CR for some EDL systems, from 
molecular dynamics (b,f) and Monte Carlo (c) simula- 
tions, are plotted in the charge reversal diagram of Fig. [5] 
We also plotted the onset of oscillations of the ionic pair 
correlation function predicted by the dressed ion theory 
(DIT) (g,h) (See table I). It can be seen that our predic- 
tion (ko« 1.260, solid line) well agrees with the onset of 
oscillations of the pair correlation function predicted by 
DIT [36[ ( fta« 1.293), for a particular case of a mono- 
valent RPM electrolyte (g). Our prediction is also in 
close agreement with the onset points of CR predicted 
by molecular dynamics (b,f) and Monte Carlo (c) simu- 
lations for several EDL systems of divalent electrolytes. 
For divalent electrolytes, our prediction is in better agree- 
ment with simulation data than DIT (h). 

In Fig. [3] u is plotted, as a function of the dimen- 
sionlcss distance X = 2x/q, for £ = 6.72, 7 = 2.4 and 
77 = 0.01,0.060.20,0.30,0.40. Note that the number and 
amplitude of oscillations increase as 77 increases. More 
importantly, however, is the fact that the absolute max- 
imum of (j(X), a*, increases by increasing 77. A similar 
behavior is observed (not shown) if 77 is constant and £ 
is increased. As we mentioned above, in Fig. [31 a* is the 
closest maximum to the wall. 

CR is not only governed by the ion-wall direct inter- 
action, but above all, it is a many body effect in which 
long range and short range correlations participate. In 
Fig. [U a* is plotted as a function of 



higher values of either parameter, £ or 77. Interestingly, 
for < 7 < 1 all curves converge to the solid line. The 
dashed and dotted lines are for 7 = 5 and 7 = 8, re- 
spectively. Curiously, the three CR curves converge for 
£ < 1.7. In general, as 7 increases the CR curves are 
shifted to the left. It should be mentioned that our dia- 
gram covers a broad range of the 77-^-7 space, including 
many cases of practical interest. 

A least squares fitting in a log-log scale was carried out 
for the data composing the solid line of Fig. O Surpris- 
ingly, we found that they satisfy the simple relationship 
[46l | = 0.0662, or equivalently 24£t7 = n 2 a 2 ~ 1.6 
, with k 2 — 2i=i Si Pii the inverse Debye length. 
This curve splits the plane in two regimes: Regime 



7 2-Ko^a 2 



(26) 



For the solid line the particles volume fraction 77 is raised 
and the Coulomb coupling parameters are kept constant 
at £ = 6.72 and 7 = 2.4. For the dashed line 77 = 0.06, 
7 = 2.4, and £ increases. Surprisingly, in both cases a* 
increases monotonically with £ in a qualitatively similar 
way. In fact, for low values of £ both curves are over- 
lapped. In this sense we refer to electrostatic and size 
correlations as symmetric. This behavior, a priori, was 
not expected since 77 and £ quantify contributions of dif- 
ferent nature, i.e., size correlations (77) and electrostatic 
correlations (£), respectively. The symmetry is also ex- 
hibited in the charge reversal curves of Fig. [21 and par- 
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FIG. 3: cr(X) for successively increasing values of r\ and con- 
stant values of £ = 6.7, and 7 = 2.4. 
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FIG. 



4: <7* as a function of £ = £77/7 = zeri/2ivaoa 2 , for 
7 = 2.4. The continuous line correspond to £ = 6.7 and 
increasing values of rj, whereas the dashed line is for r\ — 0.06 
and increasing values of £. 



ticulary for 7<1 (£r?=0.0662): when rj (£) is low, a higher 
value of £ (ry) is needed to produce charge reversal. 

From a fundamental point of view, the extra adsorp- 
tion of particles to the wall, i.e., CR, implies a gain of 
entropy respect to the bulk [1, Q, thus, adsorption is fa- 
vored. A simple way of understanding this, is through an 
analogy with the hard sphere model fluid next to a hard 
wall. In that case, size correlations induce adsorption 
of particles to the wall, which increases as the particle's 
size or concentration increase, i.e., as 77 increases. In the 
HNC/MSA theory for a RPM electrolyte, electrostatic 
and size correlations are taken into account. Perhaps 
their symmetry (Figs. [2] and 0]), can be explained by an 
effective excluded volume that mutual electrostatic re- 
pulsion between like charged ions confers them. There- 
fore, the higher value of £, the higher the particles repul- 
sion and the higher effective excluded volume, hence, the 
higher CR. 



FIG. 5: a* as a function of 7 = ^-qiCToa, with £77 = 0.220, 
in all cases. An increasing of 7 is equivalent to an increasing 
of ao- 



In Fig. [5l we evaluate the effect of 7 on <r*. For each 
curve different values of 77 and £ are combined, but a fixed 
value of £7/ = 0.220 is set. The highest considered value of 
the particles volume fraction is 77 = 0.216 with £ = 1.02. 
For this particular case, the number of adsorbed countc- 
rions per unit area, n cou , increases at a lower rate than 
00; hence a decreasing behavior of a* is observed. For 
£ = 2.04, a less pronounced, but also decreasing behav- 
ior of a* is seen. The former decreasing trend of a* vs 
7 is reversed for £ = 4.08 and more dramatically for 
£ = 16.32. Curiously, the four shown curves cross at 
7f»1.75, and below this point, a*(r/ = 1.02) <<r*(r] = 
2.04) <a*(r) = 4.08) < a*(r) = 16.32). This in turn 
implies that a* is dominated by size correlations within 
this interval. Above this interval electrostatic correla- 
tions dominate. From this perspective, electrostatic and 
size correlations do not appear to be symmetric but, by 
contrary, are unsymmctric. Both trends, increasing and 
decreasing, of the a* vs 7 curves are supported by com- 
puter simulations in studies of CR in cylindrical [4J and 
spherical geometries [l2| . 



IV. CONCLUSIONS 

We studied the electrical double layer (EDL) of a pla- 
nar wall next to a restricted primitive model electrolyte 
through a well established integral equations theory for 
inhomogeneous fluids, i.e., the hypernetted chain/mean 
spherical approximation (HNC/MSA). The HNC/MSA 
equations were derived in terms of the fundamental pa- 
rameters of the system: the particles volume fraction, 
7 ?=fPTa 3 , the ion-ion, S,=^-, and ion-surface Coulomb 
coupling, 'y=2irq/3aoa/e. To the best of our knowledge, 
for the first time a charge reversal diagram for the planar 
EDL is formally constructed. In this diagram two regimes 
of charge reversal are identified: In regime I (oscillatory) 
, charge reversal occurs for any value of 7 different from 



7 



zero and independently of the external field geometry 
whereas, in regime II (non-oscillatory), a sufficiently high 
value of 7 is required. Using numerical analysis we found 
a simple formula which distinguishes these two regimes: 
regime I is defined by the condition £77 > 0.0662, whereas 
the regime II is defined by £77 < 0.0662. This simple 
formula for estimating when charge reversal is expected 
is relevant for colloid science, since oscillations of the 
charge profiles have been associated with an effective at- 
tractive interaction between like charged macroparticles 
0, HI, Hil . The onset curve of charge reversal for the 
oscillatory regimes is in good agreement with the onset 



points of charge reversal obtained by molecular simula- 
tion for several EDL systems. Our study pointed out 
a symmetry in charge reversal between short range and 
electrostatic ionic correlations. The results of our study 
cover a broad range of practical interest. 
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